Domain Number Distribution in the Nonequilibrium Ising Model 
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We study domain distributions in the one-dimensional Ising model subject to zero-temperature 
Glauber and Kawasaki dynamics. The survival probability of a domain, S(t) ~ t~^ , and an unre- 
acted domain, Qi(t) ~ t~ s , are characterized by two independent nontrivial exponents. We develop 
an independent interval approximation that provides close estimates for many characteristics of the 
domain length and number distributions including the scaling exponents. 
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I. INTRODUCTION 

The theory of phase ordering kinetics, or domain coars- 
ening, has undergone a rapid development in recent years 
PJ. It has been established that systems quenched from 
a homogeneous high-temperature disordered state to a 
low-temperature multi-phase state do not order instanta- 
neously; instead, domains of equilibrium ordered phases 
form and grow with time as the system approaches lo- 
cal equilibrium on larger and larger scales. Generally, 
a scale-invariant morphology is developed at late times, 
and the network of domains is (statistically) independent 
of time when lengths are rescaled by a single characteris- 
tic length scale L(t), the typical domain size. This length 
scale exhibits an algebraic growth with time, L(t) ~ t v . 
However, it was recently realized that additional scal- 
ing laws characterized by nontrivial scaling exponents 
exist in such systems. Examples for such decay modes 
are the autocorrelation function, A(L) ~ L~ x ||, and 
the fraction of the system still frozen in its initial state, 
P (t) ~ t-° ||§. The latter "persistence" probability 
has since been investigated theoretically [pp|-pT| and ex- 
perimentally |H] in spin systems, interacting particles 
systems |fl p^| - |16|, Lotka-Volterra models Jn],[l8| , breath 
figures growth |19| , foams pc| , and even simple diffusion 

Similar to the domain growth exponent, v, these ad- 
ditional exponents are sensitive to the nonequilibrium 
dynamics followed by the system, and thus are funda- 
mentally different from their equilibrium counterparts. 
Precisely how many independent hidden exponents does 
a coarsening system possesses remains an open ques- 
tion. In this study, we establish that at least in one- 
dimension, additional exponents describe the survival 
probability and other more subtle statistical properties 
of domains. We examine systems with short-range inter- 
actions described by a scalar order parameter, namely the 
ID T=0 Ising model |2j| evolving according to noncon- 
served Glauber dynamics [2J] and conserved Kawasaki 
dynamics p5{ . 

This paper is organized as follows. We first define the 
domain number distribution in Sec. II. In the following 
section, we review our results for Glauber spin-flip dy- 
namics where we develop and solve analytically an Inde- 
pendent Interval Approximation (IIA) that assumes no 



correlations between adjacent domains. The IIA predic- 
tions compare well with Monte Carlo simulations by giv- 
ing a correct description of the domain statistics as well as 
good estimates for the underlying exponents. In Sec. IV 
we show that nontrivial exponents underly the zero tem- 
perature limit of the ID Ising model with Kawasaki spin- 
exchange dynamics as well. The IIA, when carefully 
modified to conserved dynamics, turns out to be equally 
useful in this case. Summary and conclusions are given 
in Sec. IV. 



II. DOMAIN NUMBER DISTRIBUTION 

Although we focus in this study on the Ising model, the 
statistical properties of domains we are concerned with 
are relevant to arbitrary coarsening processes in one spa- 
tial dimension. For example, we ask, what is the domain 
survival probability S(t), i.e., the probability that a do- 
main, initially present at the system at time t = 0, is still 
present at time t (see Fig. 1). We will present theoreti- 
cal and numerical evidence supporting an algebraic long 
time decay of this survival probability, 



(1) 



Such a behavior is robust, as the exponent ip is not sen- 
sitive to the initial state of the system (provided long 
ranged correlations are absent). Our results will also 
strongly suggest that the exponent tp is nontrivial, i.e., 
it cannot be extracted from so-far known exponents as- 
sociated with the Ising model. 

+ - + - + + 



m=3 / m=\\ 

Fig. 1. Domain motion in the Ising-Glauber model. Sur- 
viving domains are marked by +, annihilated domains by — . 
The domain number at a later time is also indicated. 
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In principle, a surviving domain may undergo coales- 
cence with other similar phase domains. Thus, a natu- 
ral generalization of the domain survival probability is 
Qm(t), the density of domains composed of m original 
domains (see Fig. 1). This quantity satisfies the initial 
condition Q m (0) = <5 TOj i. The total domain density, N(t), 
is given by N(t) — ^ m Q m (t), while the domain survival 
probability counts initial domains that have not shrunk 
and hence contains the density Q m (t) with weight m 

S(t) =J2mQ m (t). (2) 

m 

The average number of domains contained within a sur- 
viving domain (m(t)) = S(t)/N(t) grows algebraically 
according to (m(t)) ~ t v ~^ with v the domain decay ex- 
ponent, N(t) ~ t~ v . If the behavior of Q m (t) is truly 
self-similar, it should follow the scaling form 

Q m (t)~t^ 2 »Q(mt^). (3) 

The scaling function Q(z) exhibits the following ex- 
tremal behavior 

Q W ~ { expi-Kz") z>L ^ 

The small argument tail describes domains that contain 
a very small number of initial domains. In particular, the 
quantity 

Qx{t)~r s (5) 

is of special interest: It gives the density of domains 
which avoided merging with their neighboring domains 
up to time t. 

The inequalities Q\(t) < ^ m Q m («) < 2 m roQ m (t) 
lead to the bounds ip < v < S. Taking into account 
that at least one surviving domain surrounds a persis- 
tent spin gives P(t) < S(t), where P(t) ~ t~ e is the 
density of persistent spins. Thus we arrive at another 
upper bound if> < 9 for the exponent ip. These bounds 
suggest that the domain decay rate is the slowest in the 
problem. We shall show below that these bounds are 
strict for the Ising model and more generally for the q- 
state Potts model. However, for the Potts model with 
q — > 1 or q — > co, and for a few other models p6|] some 
of these exponents are equal to each other. 

A useful relation between the scaling exponents can be 
obtained by substituting m = 1 in Eq. (^) 

8-v=(v-^){l + a). (6) 

Thus, among the three exponents ip, <5, and er, only two 
are independent. It is well known that under noncon- 
served (conserved) dynamics i/=l/2(i^ = l/3) [0. 

Quite obviously, domains disappear when their size 
vanishes, and therefore the domain size and number dis- 
tributions are intimately related. Thus, domain survival 
properties involve the distribution of domains of size 



n consisting of m original domains at time t, denoted 
by P n ,m{t)- The aforementioned number distribution is 
Qm(t) = ^ n P n ,m{t), and consequently, the domain sur- 
vival probability is S(t) = J2 n .m rnP ",m(t)- 

As will be seen later, studying the joint size-number 
distribution requires detailed knowledge of the domain 
size distribution P n (t) — ^ m P n ,m(t). This distribution 
obeys the normalization conditions 

l = 5>P n (t), N(t)=J2Pn(t). (7) 

n n 

Length conservation implies the first relation, while the 
second relation gives the total domain density. Since the 
average domain length grows as n ~ t" , the length dis- 
tribution follows the scaling form 

P n (t)~t- 2 »V(nt-")- (8) 

All of the above scaling behavior emerges from the ap- 
proximation detailed below. Furthermore, it is satisfied 
by the simulation data. In the next section, we develop 
an approximation scheme that helps elucidate many of 
the qualitative and quantitative features of the domain 
size and number distributions. 



III. NONCONSERVED GLAUBER DYNAMICS 



We start with the ID Ising model subject to T = 
Glauber dynamics 0] . To examine the role of the num- 
ber of equilibrium phases we also consider a generaliza- 
tion of the Ising model, the g-state Potts model. In higher 
dimensions, the g-state Potts model is relevant to phys- 
ical situations for q = 2 (the Ising model) and addition- 
ally for q = 3,4, oo ^?J. For instance, the q — oo case 
describes several cellular structures ]28| including poly- 
crystals p3] , foams [^0| , soap froth |30fl , and magnetic 
bubbles |31fl . 

We consider uncorrelated initial conditions where each 
of the q phases is present with equal density l/q. The 
T = Glauber-Potts dynamics proceeds by selecting a 
spin at random and changing its value to that of one of 
its randomly selected neighbors. Thus, domain walls per- 
form a random walk and upon contact, they annihilate 
or coalesce, depending on the state of the corresponding 
domains [p2|-p4[. Identifying a domain wall with a par- 
ticle, (A), and absence of a domain wall with a hole (0), 
one finds the single-species diffusion-reaction process 

1 1 q-2 

AQ^OA, AA^\ 00, AA ^> AO or OA. (9) 

The rates indicate the relative probabilities by which each 
event occurs. 
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A. Domain Size Distribution 

Ignoring correlations between neighboring domains al- 
lows us to develop an approximate theory for the time- 
evolution of the domain distribution. Approximations 
that are similar in nature proved useful in studies of re- 
lated reaction-diffusion processes J2l],|32],[55| . 

The joint number distribution requires knowledge of 
the length distribution and we start by deriving a mas- 
ter equation for P n (t). Under the assumption that the 
lengths of neighboring intervals are uncorrelated, we 
write the following rate equation [Eel 



dP n 
dt 



P-n-l + P-n+l ~ %Pn 
~n-2 

^PiPn-^i-NiPn + Pn-i) 



(10) 



Pi 



(q- l)N 2 



with N{t) = P n (t) the total domain density and the 
boundary condition Po(t) = 0. The first three terms 
reflect that domain walls perform a random walk with 
hopping rate set to 1/2 without loss of generality. The 
last two terms are due to domain annihilation: the con- 
volution term accounts for domain merger and the last 
term for domain loss. In the g-state Potts model collision 
of domain walls results in annihilation with probability 
or in coalescence with probability ^E§- O nrv anni- 
hilation events affect the domain distribution and thus 
the prefactor of the annihilation terms. Using the 
sum rules of Eq. @, one verifies that the total length is 
conserved and the total domain density decays according 
to the exact rate equation 



dN 
~dt 



Pi- 



(11) 



The diffusion term in Eq. (0) implies (n(t)) ~ t 1 / 2 , 
and since (n) ~ A -1 the correct decay exponent v = 1/2 
p4| is recovered. In the following, we will need to de- 
termine the asymptotic prefactor A, N(t) ~ At^ 1 ^ 2 , 
A = f dxV(x), with the scaling function V(x) defined 
according to Eq. (|||). The density rate equation (0) im- 
plies Pi ~ V'{0)t- 3 ^ 2 with V'(Q) = ^-A. 

A quantitative analysis of Eq. (|l(]) may be carried by 
treating the variable n as continuous. The quantity V(x) 
satisfies 



where V = dV/dx and V * V = f* dyV{y)V{x - y). 
The normalized Laplace transform of the scaling func- 
tion V(x), p(s) = A~ x J °° dx e~ sx V(x), obeys 



dp 
ds 



El 

qs 



2s 



q-2\ q 
P 

qs 



1 



qs 



(13) 



subject to the boundary condition p(0) = 1. The trans- 
formation's) = l~qs 2 — qs-j^ lny(s) reduces the Riccati 
equation (O) into the parabolic cylinder equation, 



ds 2 



2 

1 + - 

q 



y = o. 



(14) 



The solution to Q reads y(s) = C-D 1/q (-sV2) + 
C + D 1 / 9 (sv / 2) J with D 1 / q (x) the parabolic cylinder func- 
tion of order 1/q p&. The large s behavior of p(s), 



p(s) ~ ^2q~ s 2 ' implies C_ = 0, and we get 



P(s) = l-qs 2 



qs—bxDi/ q (sV2). 



(15) 



The normalization condition nP n (t) = 1 can be re- 
duced to Ap'(Q) = — 1. This allows us to determine the 
constant 



.4 



r[i-£] 



X 12 2q\ 



(16) 



where T denotes the gamma function. In deriving (|l6| ) 
we have used the properties |3§] 



D c {x) ~ ;r c exp(-a; 2 /4)[l + C , (a;" 2 )], 



(17) 



and 



" D'M = - J ■ ( 18 ) 



r(l/2-c/2)' 



r(- c /2) 



The value of the constant A predicted by the IIA may 
be compared to the exact one, A oxact = (1 — q^ 1 )/^ |3J] 
In the extreme cases of q = 1 and q = oo the prefactor 
A is exact. The mismatch is worst (roughly 20%) for the 
Ising (q = 2) case where A = r(3/4]/T(l/4) = 0.337989 
while A cxact = (47T)- 1 / 2 = 0.28209 |J. 

The IIA predicts the correct qualitative behavior of 
length distribution in the limits of small and large inter- 
vals 



V(x) 



2q 



X < 1, 



qAexp(~Xx) x ^> 1. 



(19) 



The linear small size behavior is seen from the large s be- 



havior p(s) 



(g-i) 

2q 



s 2 . On the other hand, the exponen- 



tial tail follows from the behavior of the Laplace trans- 
form p(s) ~ q\j (s + A) near its pole at negative s = — A, 
given by the first zero of D 1 / q (— A\/2) = 0. For the Ising 
case one has A = 0.5409. This value should be com- 
pared with the exact value A = ((3/2)/ 4^ = 0.368468 
obtained by Derrida and Zeitak |37j and the approxi- 
mate value A = 0.35783 obtained by Alemany and ben- 
Avraham |35|. 
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B. Domain Size-Number Distribution 

We are now in a position to tackle the joint size-number 
distribution, P n . m (t), which captures both the spatial 
and "historical" characteristics of the coarsening domain 
mosaic. The corresponding rate equation is a generaliza- 
tion of Eq. © 



dP 
dt 



= Pn- 



Pn 



+ l,m 



-2P„ 



(20) 



Pi 



(q~l)N 2 



^ ^ Pj,jPn— 1 — i,m— j ^{,Pn.m ~\~ Pri—X,m) 



with the initial condition P n . m (0) = <5 n ,i<5m.i and the 
boundary condition Po,m(i) = 0. The variable m is al- 
most mute as it appears in a nontrivial way only in the 
convolution term. One should verify that this master 
equation is self-consistent. First, by summing over to, 
we recover Eq. (|l0|). Second, it implies that the do- 
main survival probability satisfies the exact linear equa- 
tion dS/dt = - J2 m m Pl,m- 

We have not succeeded in solving for the joint distri- 
bution. Nevertheless, it is possible to obtain analytically 
many interesting properties of Eqs. (p0|), including the 
scaling exponents. Given Eqs. (|20| ) are recursive in to, 
one can try to solve for P n> i(i), then for P n ,2{t), etc. A 
solution for the former quantity already allows to deter- 
mine the scaling exponent 5. Thus let us consider the dis- 
tribution of domains which have not merged with other 
domains up to t, R n (t) = P n ,i(t). For such domains, the 
convolution term vanishes and they evolve according to 
the linear rate equation 



dR n 
~dt 



Rn-l + Rn+1 — 2i?„ 



Pi 



(q-l)N 



(R n + Rn 



(21) 



with the initial condition P„(0) = 8 nt i and the bound- 
ary condition Ro(t) — 0. In the continuum limit we 
again replace P„_i + R n +i — 2i?„ by d 2 R/dn 2 and 
R n + R n -i by 2R n to find a diffusion-convection equa- 
tion for R n (t). The transformation R n — + R n N 2 / q re- 
duces this equation to the diffusion equation for R n , 
which is solved to yield R n (t) ~ 7V 2 /«t~ 1 ^(nt" 1 / 2 ), with 
IZ(x) = xexp(— x 2 /4)/y / 7r. The total density of unre- 
actcd domains is Qi(t) — ^2 n Rn ~ t _3 ", which gives 
the decay exponent 



r 1 1 
S =2 + ~ q 



(22) 



Obtaining the second independent exponent -0 is more 
involved. The natural approach, i.e., a direct inves- 
tigation of the domain number distribution Q mi ap- 
pears to be useless, as it requires knowledge of Pi ym and 
hence the entire P n ,m- The domain survival probabil- 
ity can be alternatively obtained by considering U n (t) = 
Em m ^«^W- Tnis quantity obeys 



dt 



Un~l + Un+1 — 2C/„ 
n—2 



(23) 



Pi 



(q - 1)N 2 



2j2U t P n -l-r ~ N(U n + U n -l) 



i=l 



obtained by summing Eqs. fl20|). We write U n (t) in a 
scaling form U n (t) ~ i"^" 1 / 2 ^™^ 1 / 2 ). Asymptotically, 
the domain survival probability reads Sit) ~ Bt~^ with 
B = J dxU(x). The scaling distribution satisfies 



U" + \{xU)'+(^--\ 



u 



1 



U*V = 0. 



(24) 



The normalized Laplace transform of the scaling function 
U{x), u(s) = B^ 1 / °° dx e~ sx U(x), obeys 



du ( p(s) + qjj - 1 , \ 20 

— = 2 \-s \u , 

ds \ qs ) s 



(25) 



and u(0) = 1. In deriving (25) we used the relation 
U'(0) = Bip, found by integration of Eq. (pjj), combined 
with A = J dxV(x). Substituting the explicit expression 
( |l5| ) for p(s) into Eq. (|25|), and solving for u(s) yields 

/oo 
drr-^- l Dl /q {ry/2). (26) 

This solution is consistent with the anticipated s — » oo 
behavior, u(s) ~ i/js~~ 2 . Furthermore, evaluating Eq. ( p6| ) 
near the origin gives u(s) = 1 + F(ip)s 2 ^ + Cs + ■ ■ ■. 
Therefore, for u'(s) to be finite near s = 0, we must have 
F(ip) = 0. Evaluating F{ip) gives 



0= / drr-^D 1/q (r)D[ /q (r), 



(27) 



an eigenvalue problem that can be solved numerically to 
obtain the exponent ip (see Table 1). In the most in- 
teresting case of integer q the domain decay exponent ip 
appears to be irrational, in contrast with S. 

It is useful to consider the limiting cases that turn 
out to be solvable. The q = oo limit is especially sim- 
ple |3^| as only domain walls coalesce but cannot an- 
nihilate and thereore Eq. (BQ) is linear and thus exact. 
Furthermore, the domain size number distribution fac- 
torizes, P n>m — P n (t)5 m ^i, since similar phase domains 
never coalesce and therefore the domain number is triv- 
ial, to = 1. Thus N(t) = S(t) = Qx{t) ~ (Trf)- 1 / 2 and 
v = ip = S = 1/2. Additionally, the scaling function is 
P(x) = xexp(-x 2 /4)/V^. 

Before going to the opposite limit q — > 1, we first note 
that the Potts model with arbitrary q > 1 can be mapped 
onto the Ising model with magnetization /i = 2/q — I. 
Thus the q — > 1 limit corresponds to the vanishing vol- 
ume fraction of minority domains. Therefore minority 
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domains cannot "meet" , so majority domains change ap- 
preciably only due to coalescence. Thus they never dis- 
appear, i.e., S(t) — 1 and ip = 0. A majority domain 
remains unreacted if both its minority neighbors sur- 
vive, implying Q\(t) — N 2 (t) and 8 = 2v = 1. In 
the above rate equation description, the diffusion term 
becomes negligible and the II A is exact. The scal- 
ing functions V (x) and Q(z) are identical exponential 
functions. However, unlike to the q = oo case, the 
joint distribution is not a product of the single func- 
tion variables pq . The joint size- number distribution 
still obeys the scaling law, P n ,m(t) ~ t~ 5 ' 4 $(x, y), with 



exp(— x — y/2x). The scaling variables 



<J>(a;, y) = x^ 1 / 2 
x and y are quite unusual, though [|38[|26 |. Indeed, in- 
stead of the naive scaling variables mV^ 2 and ni -1 / 2 , 
one has x = (m + n)(7rf)~ 1//2 and y = (m — n)(7rf)~ 1//4 
p6fl . The former scaling variable x is just the sum of the 
naive scaling variables, while the latter "diffusive" scale 
y is hidden. This suggests that generally for the g-state 
Potts model with q < oo the joint distribution is not nec- 
essarily the product of single variable functions and the 
scaling, if holds, may be rather different from the naive 
form with scaling variables nt~ v and mt^'~ lJ . 



C. Simulation Results 

To test the IIA predictions, we performed numerical 
simulations on a spin chain of size L = 10 7 . Random 
initial conditions and periodic boundary conditions were 
used. The simulation data represents an average over 10 
different realizations. For the Ising case, we found the 
exponent values ip — 0.126(1) and 8 — 1.27(2) (see Fig. 
2). These values should be compared with the IIA pre- 
dictions of ip = 0.136612 and 8=1. The IIA neglects 
correlations that do build up between neighboring do- 
mains and thus is not exact. Furthermore, the effects of 
the correlations is nontrivial, as one exponent is smaller 
than predicted while the other is larger. As was the case 
for the persistence exponent, 9, the domain exponents 
strongly depend on q. Numerical values of the exponents 
ip and 8 are summarized in Table 1 for representative val- 
ues of q. As q increases, the approximation improves and 
eventually becomes exact for the extreme case q = oo. 
Thus, ip is overestimated by up to 10% and <5 is underes- 
timated by up to 25%. 

We performed several checks to verify that the asymp- 
totic behaviors of Eqs. (|Tj) and (||) are robust. For exam- 
ple, they are independent of the initial domain wall con- 
centration (provided that the correlations in the initial 
condition are short range). We conclude that ip and 8 are 
nontrivial exponent, i.e., they cannot be extracted from 
the known exponents associated with the Ising-Glauber 
model. Similar to the persistence exponent, 0(g), the ex- 
ponents appear to be irrational except for the limiting 
cases q = oo {ip = S = |, a = 0) and, maybe, for q = 2 
(i> = l,5=l<r = l). 



The numerical simulations also confirm that the distri- 
bution function Q m (t) scales according to Eq. (||). The 
scaling function Q{z), defined in Eq. (||), decays expo- 
nentially for large argument (p = 1) and is algebraic 
for small argument. The scaling relations combined with 
the simulation values give a — 1.05(5). This is consis- 
tent with the linear behavior seen in our simulations for 
z«l. Comparing with Eq. ([p]), we conclude that sim- 
ilar scaling functions underlie the domain number and 
size distributions in the q — 2 case. 
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Fig. 2. Monte Carlo data for the Ising-Glauber model. The 
domain survival probability S(t), the domain density N(t), 
and the density of unreacted domains Qi(t) are shown (top 
to bottom). The inset plots the local slope — din S(t)/d Int. 
Typically, it is stable over a large temporal range, and thus 
can be used to find the scaling exponents and to estimate the 
error, typically of the order 0.001. 

On the other hand, direct numerical integration of 
Eq. ( |20| ) reveals a number distribution, Q m (t) , that scales 
according to Eq. (|^) , and has an exponential tail in agree- 
ment with the simulation results. Moreover, the emerg- 
ing S(t) falls within 5% of the actual survival probability 
over a significant temporal range, t < 10 3 . In summary, 
in addition to predicting the correct scaling behavior, 
Eq. (po| ) provides a good approximation for many quan- 
titative features of the domain distribution, and in par- 
ticular, good estimates for the decay exponents. 



MC 



<1 


i> 


8 


(T 


2 


0.126 


1.27 


1.05 


3 


0.213 


0.98 


0.67 


4 


0.267 


0.85 


0.50 


8 


0.367 


0.665 


0.24 


50 


0.476 


0.525 


0.03 


00 


1/2 


1/2 






Eq. (20) 



0.136612 
0.231139 
0.287602 
0.385019 
0.480274 
1/2 



1 

5/6 
3/4 
5/8 
13/25 
1/2 



Table 1: Domain exponents for the g-state Potts model in 
one dimension. Local slopes analysis was applied to the simu- 
lation data. The theoretical ip is from Eq. 



and 5 — | 
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IV. CONSERVED KAWASAKI DYNAMICS 

We turn now to applying the above methods to the con- 
served counterpart, the spin-exchange Kawasaki dynam- 
ics , which describe spinodal decomposition in binary 
alloys and phase separation in binary liquids. Although 
some qualitative features are known p9|-^l|, theoretical 
understanding of the Ising-Kawasaki model is still far 
from complete even in one dimension. In the following 
we limit ourselves to the two-phase Ising case. 



A. Reduction to Domain Diffusion 



Thus, the appropriate zero-temperature limit of the 
Ising-Kawasaki model is realized by taking the limits of 
infinite "physical" time t p h ys — ► oo, while keeping the 
modified time t = i p hy S exp[— 4 J/T] finite. Hence, en- 
tire domains perform a random walk with rate inversely 
proportional to their lengths (we ignore an anomaly con- 
cerning domains of length 2 as it is irrelevant asymptoti- 
cally). Heuristically, it may be argued that as diffusion is 
the primary coarsening mechanism, the following scaling 
for the average domain size, L ~ yDt, holds. However, 
since the diffusion coefficient and the domain size are re- 
ciprocal, D ~ we obtain L ~ i 1 / 3 in agreement with 
the well-known behavior of systems with conserved scalar 
order parameter (l). 



We start by formulating the appropriate zero- 
temperature limit of the Ising-Kawasaki model. Con- 
sider a two- phase system, e.g., the Ising model (spins 
up and down) or a binary alloy (atoms of type A 
and B). At zero temperature energy raising transi- 
tions are forbidden and only two moves are allowed: 
the energy-decreasing "coarsening" transitions ABAB — > 
AABB and the energy-conserving "diffusion" transitions 
ABAA — > AABA. This dynamics ultimately drives the 
system to a frozen configuration consisting of strings of 
alternating domains each of length > 2 which could not 
evolve further Ji2|| . This "jamming" behavior arises from 
the nonergodic nature of the zero-temperature Kawasaki 
dynamics and it is a robust one: It is independent of 
the relative transition rates , as well as the spatial 

dimension 

Thus, to recover sensible coarsening, one must con- 
sider the zero-temperature limit. Let us assume that 
temperature is positive; when it is sufficiently low, T <C J 
where J is the exchange coupling, the correlation length 
£ ~ e J / T is very large and therefore the system exhibits 
coarsening as long as the mean domain size is small com- 
pared with the correlation length. Below, we focus on 
this intermediate-time regime where the description is 
drastically simplified j3^,Q. We assume that the initial 
stage has been already completed so that single-spin do- 
mains disappeared. Coarsening will occur only when a 
spin splits off a domain wall and penetrates a neighbor- 
ing domain (say of size L). The splitting process occurs 
with a very small rate exp[— A J/T]. Then, this spin dif- 
fuses inside the domain until it is eventually adsorbed by 
its boundaries. The corresponding probabilities are well- 
known from elementary probability theory p6| . The spin 
will be absorbed by the boundary from which it was is- 
sued with probability 1 — 1/ L. This spin may also be ab- 
sorbed by the opposite boundary resulting in a one lattice 
site hop of the entire domain. Thus, the hopping rate is 
L" 1 exp[— 4J/T]. Rescaling time, t — ► £exp[— 4J/T], the 
spin diffusion will proceed with a huge rate exp[4J/T] 
and it therefore may be treated as instantaneous while 
domain hops proceed with a finite rate reciprocal to the 
domain size. 



B. Domain Size Distribution 

For simplicity, we consider the case where the two equi- 
librium phases are equivalent, as is the case for random 
initial conditions. Modifying Eq. ( |l0|) to account for do- 
main diffusion, and assuming neighboring domains are 
uncorrelated, the domain size distribution evolves ac- 
cording to 



dPn 

dt 



— L (P n —i — 2P n + P n +\) 



(28) 



Pi_ 

N 2 



r X PiPj ~ Nin- 1 + L-^Pn 

i+j=n 



with the domain density N and the inverse average do- 
main size L~ x defined via 



N = f^P n , L-i = (n-i) = ^ n ~ 1Pn . 

„. = 1 l^n 



(29) 



The diffusion term in Eq. (|28|) accounts for change in n 
due to hoping of a neighboring interval. The convolu- 
tion term accounts for gain due to domain merger, and 
the last two terms represent loss due to domain colli- 
sion or domain merger. This rate equation conserves 
the total length, ^ n nP n = 1, and summation over n 
shows that the total domain density evolves according to 
N = -2£- 1 Pi. 

The length distribution scales according to Eq. (|§), 
P n (t) ~ t~ 2v V{nt~ v \ with the correct scaling exponent 
v = 1/3. By inserting that scaling form into Eq. ( Ba ) we 
arrive at an integro-differential equation for V{x) which 
is very cumbersome as it involves yet unknown moments 
of the distribution. We thus resort to numerical in- 
tegration of Eq. (p8|). The results compare well with 
Monte-Carlo simulations of the Ising-Kawasaki model 
(see Fig. 3). For example the estimate for the asymptotic 
prefactor A (defined via N(t) ~ At^ 1 ^ 3 ), falls within 
roughly 5% of the actual value: A M c = 0.441 ± 0.001 
while Aha = 0.415 ± 0.005. 



G 



0.40 



0.30 



0.20 



0.10 



0.00 



MC 
IIA 



0.0 



200.0 



400.0 



600.0 



800.0 



1000.0 



Fig. 3. Domain density in the Ising-Kawasaki model. The 
Monte Carlo simulation data (MC) represents an average over 
10 systems of size 10 5 . The IIA was obtained by integrating 
Eq. ( |28| ) numerically. 

The length distribution emerging from the IIA has the 
same limiting behavior as in the nonconserved case, i.e., 
it is linear at small n and exponential at large n. While 
the former agrees with our simulation results, there is a 
disagreement for the latter. Our data is consistent with 
a Gaussian tail, i.e., V{x) ~ exp(— x 2 ), for x ^> 1. 



C. Domain Size-Number Distribution 

Given the results in the conserved case, it is natural 
to study the domain exponents and to examine the use- 
fulness of the IIA approach in the conserved dynamics 
case. In analogy with Eq. (|20|), the master equation for 
the domain size-number distribution is written 
dP 



dt 
N 2 



(30) 



J2 r'-Pi.kPjj-Nin-i+L-^Pn 



i+j—n k+l—m 



Summing the above equations over m, we indeed recover 
Eq. (^||) for the length distribution. 

To determine the exponents, it is again simpler to con- 
sider the distributions R n (t) and U n (t) instead of the 
joint distribution P„ jm (i). The density of single parent 
domains, R n (t) = P n: i(t) : evolves according to the linear 
rate equation similar to Eq. (|2l| ) 



dR n 
~df 



L 1 {R n +l + Rn-l — 2i? n 

^{n^+L- l )R n . 



(31) 



We expect that the distribution R n (t) scales according to 
R n (t) ~ i _5_1 / 3 7^(nf -1 / 3 ). Integrating this equation we 



get Qi(i) = J2 R n{t) ~ t~ s with 5 = 0.645. On the other 
hand, Monte-Carlo simulations of the domain diffusion 
process give S = 0.705 (see Fig. 4). 

The domain survival probability can be found by us- 
ing the auxiliary function U n (t) — J2m m Pn,m(t) which 
satisfies the analog of Eq. ( ^3|) 



dU n 
dt 



L-\U n+1 + U n - 1 -2U n 



(32) 



P_ 

'N 2 



2 i^PiUj - Nin- 1 + L~ l )U n 

i+j—n 



This distribution U n (t) should scale according to U n (t) ~ 
t~^~ x IHA.{nt~ x l 3 ). The domain survival probability 
S(t) = Yin Un(t) then decays according to S(t) ~ t~^ . 
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Fig. 4. Monte Carlo data for the Ising-Kawasaki model from 
the same simulation as in Fig. (3) . The domain survival prob- 
ability S(t), the domain density N(t), and the density of un- 
reached domains Q\(t) are shown (top to bottom). 

Again, the agreement with the simulation is remark- 
able. Numerical integration data give an estimate of 
if> = 0.147 while Monte-Carlo simulations (see Fig. 4)) 
give ij) = 0.130. We also verified that the scaling relations 
of Eqs. (||) and (^) are satisfied by the IIA as well as the 
simulation data. We conclude that in the Kawasaki case 
as well, nontrivial exponents characterize domain statis- 
tics. Furthermore, the approximate approach reproduces 
most qualitative features of the domain size and number 
distribution, and provides good estimates for the scaling 
exponents. 

We now describe how to obtain tp in the limit where one 
of the two phases occupies a vanishing volume fraction. 
Denote by La(*) and Ls(t) the average sizes of minority 
and majority domains, respectively. They both grow as 
i 1 / 3 but remain greatly different throughout the evolu- 
tion, La <C Lb- The domains diffusion rates, Da ~ L^ 1 
and Db ~ Lg , thus greatly differ as well: Da Db- 
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In principle, two neighboring minority domains can over- 
take a separation distance of order Lb and coalesce; this 
requires the coalescence time t c ~ L b /Da ~ LaL b . 
On the other case, a minority domain can shrink due to 
diffusion of neighboring majority domains; this requires 
the shrinking time t s ~ L\/Db ~ L\Lb- We see that 
t s *C £ c; so we should just take into account the shrinking 
of minority domains. 

Thus majority domains do not disappear, implying 
S(t) = 1 and ip = 0. For the minority phase we antici- 
pate ip — v — while for the symmetric case of equal 
concentrations ip = 0.130. This indicates that similar to 
the nonconserved dynamics case, the domain exponents 
vary continuously as the volume fraction is varied. 



V. CONCLUSIONS 

In summary, we investigated the one-dimensional 
Ising model subject to zero temperature Glauber and 
Kawasaki dynamics. We introduced the domain size- 
number distribution and showed that it obeys scaling 
and is characterized by two independent nontrivial decay 
exponents. Similar to the persistence exponent, these 
exponents are sensitive to the type of the dynamics and 
the volume fraction of the (globally or locally) conserved 
equilibrium phase. We also introduced an approximation 
which is based on terminating the hierarchy of rate equa- 
tions describing the domain density. This approximation 
is very useful in predicting the qualitative nature of the 
domain distribution as well as estimating important pa- 
rameters including the scaling exponents. In the proper 
T=0 limit of the Ising model with Kawasaki dynamics, 
this approximation is especially useful as very little is 
known analytically about the domain distribution. 

It will be interesting to generalize the domain survival 
concept to higher dimensions. At least for the q — > oo 
limit of the Potts model, domains are well defined, and 
such a generalization is possible. The nonconserved dy- 
namics can indeed be studied in this limit, and the do- 
main exponents ip — 6 = S = d/2 (for d > 2) have been 
reported [ p6[ , consistent with simulations pC| and with 
experiments on d = 2 soap froths (4?j . 

Recently, it was pointed out that coarsening mosaics 
may be characterized by more than one algebraically 
growing length scale, and that morphologies consisting 
of domain and super-domains may exist [p.7| . This, to- 
gether with the above results suggest that our current 
understanding of such systems is only partial. Domain 
statistics indicates that several nontrivial decay laws un- 
derlie the evolution of elementary processes such as the 
nonequilibrium Ising model. These nontrivial exponents 
do not emerge naturally from studies of traditional quan- 
tities such as spatiotemporal correlations. It remains a 
challenge to find and obtain these underlying "hidden" 
exponents from a more systematic method. It is also in- 
triguing whether an entire hierarchy or a finite number of 



independent decay modes are present in these systems. 



This paper is dedicated to Leo Kadanoff on the occa- 
sion of his 60th birthday. The research of PLK has been 
supported by grants from NSF and ARO. 
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